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Abstract 

We  consider  a symmetric  boundary  integral  formulation  associated  with  a mixed  bound- 
ary value  problem  defined  on  a domain  Cl  € IR2  with  piecewise  smooth  boundary  F. 
We  assume  that  Cl  is  mapped  onto  itself  by  a finite  group  Q of  congruences  having  at 
least  two  distinct  elements.  Hence,  we  can  decompose  the  related  symmetric  Galerkin 
BEM  problem  into  independent  subproblems  of  reduced  dimension  with  respect  to  the 
complete  one.  Shape  functions  for  each  subproblem  can  be  obtained  from  classical  BEM 
basis,  ordered  as  a vector,  applying  suitable  restriction  matrices  constructed  starting 
from  group  representation  theory. 


1 Introduction 

Let  Cl  C 1R2,  be  a bounded  domain  with  a piecewise  smooth  boundary  T.  The  boundary 
T is  partitioned  into  two  non  intersecting  open  subset  and  T2,  with  F = l’i  (J  F2  = 
UEi  r,  P being  an  open  strainght  line  segments.  In  the  following  we  always  assume 

(1.1) 
(1.2) 

(1.3) 

In  (1.1)  !/(•)  is  an  elliptic  partial  differential  operator  of  second  order,  U(x,y)  its  fun- 
damental solution  (see  [4]  for  a general  discussion).  In  (1.2)  §)(  denotes  the  derivative 
with  respect  to  the  outher  normal  n to  T,  and  u*  and  q*  are  given  functions.  Applica- 
tions of  (1.1)-(1.2)  are,  for  instance,  boundary  value  problems  in  potential  theory  and 
in  elastostatic.  From  (1.3)  it  is  clear  that  if  we  want  to  recover  u in  $7  we  have  firstly  to 
know  the  remaining  Cauchy  data,  since  in  (1.2)  these  functions  are  given  only  partially. 
Taking  the  limit  of  u(x)  for  x € Ta  and  the  normal  derivative  |^(x)  for  x E T2  in  this 
formula  and  using  the  jump  relations,  one  finds  the  system  [2] 

f U(x,y)q(y)dy  - / ~^—U{x,y)u{y)dy  = fi{x),  leFj, 

J ra  J r2  any 


— . — * 


measTi  > 0.  The  solution  of  the  mixed  boundary  value  problem 

L(x)u(x)  — 0 in  Cl, 

3x1 

u(x)=u*(x)  on  Ti , q(x)  :=  — = q*(x)  on  T2, 
can  be  expressed  by  the  representation  formula 

u(x)  = J U(x,y)q(y)dy  - J -^-U(x,y)u{y)  dy,  xeCl. 
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- £ ^u(x,y)q(y)dy+  £ ■^^~U(x,y)u{y)dy  = f2(x),  xeT2.  (1.4) 

In  order  to  perform  the  Galerkin  method,  we  need  a family  of  finite-dimensional  sub- 
spaces {Uh,p{ r)}  defined  on  F.  Let  us  define  a mesh  for  each  T-3:  = (J^i  F^ 

such  that  r3  f is  an  open  segment.  We  define  for  p > 0,  h > 0,  Uh,p{T\)  to  be  the  set 
of  functions  on  Fi  whose  restrictions  to  C Fi  belong  to  the  set  of  all  polynomials  of 
degree  < p on  T3h  i.  Moreover,  for  p > 1,  U^  p(T2)  will  denote  those  continuous  functions 
on  T2  whose  restrictions  to  C r2  belong  to  C'0(r2)  and  which  vanish  at  the  end  points 
of  r2.  The  approximating  boundary  element  shape  functions  of  degree  p > 0 are  defined 
through  the  standard  assembling  of  the  local  basis  functions  defined  on  each  F)t  ..  We 
then  define 

Uh,P(r ) '■=  sPan  ipt)  : <fii  € U^p(T2),  ipe  6 UhtP{Ti)}.  (1.5) 

The  corresponding  symmetric  Galerkin  boundary  elements  scheme  for  (1.4)  leads  to  a 
linear  system  of  the  form 

A£=b.  (1.6) 

If  the  boundary  F presents  symmetry  properties,  we  will  exploit  them  to  reduce  the 
computational  cost  of  the  solution  of  (1.6),  using  a decomposition  result  for  the  Galerkin 
boundary  element  problem  that  we  will  introduce  at  the  end  of  the  next  section. 

2 Matrix  representation  of  a finite  group  of  congruences  and 
projection  operators 

Let  Q be  a finite  group  of  t congruences  ( t > 2)  of  the  Euclidean  space  IT"  (m  = 2,  3). 
The  group  Q can  be  described  by  orthogonal  matrices  7,  of  order  m.  Let  {71, . . . , 74}  be 
the  elements  of  G,  7i  the  identity  matrix.  From  the  theory  of  group  representation  [5]  it 
follows  that  any  finite  group  Q admits  a finite  number  q of  unitary  irreducible,  pairwise 
inequivalent  matrix  representations 

{w(1)(7i)>,  {w(2)(7i)},...,  {w(?)(7i)}  (*  = 1, ...,«) . (2.1) 

Let  dp  be  the  order  of  the  representation  (0/^(7,)},  i.e.,  the  order  of  the  matrices 
(7 »)•  The  number  q of  the  representations  (2.1)  and  the  orders  d\, . . . , dq  only  de- 
pend on  Q.  Any  representation  {0/^(7*)}  of  order  de  > 2,  can  be  replaced,  in  the 
system  (2.1),  by  an  equivalent  unitary  representation.  Representations  of  order  1 are 
univocally  determined.  We  observe  that,  if  7 , and  7 are  two  elements  of  Q,  then 
w(^)(7i7j)  = (7j)  > = [£*>^(74)]*,  where  [o>^(7,)]*  denote  the  trans- 

pose of  the  matrix  0/^(7*).  Always  from  the  theory  of  group  representation  it  follows 

that  q < t and  the  relation  d\+  d\- 1 + d2q=t  holds.  Furthermore,  q = t if  and  only 

if  di  = d2  = • • • = dq  — 1.  Having  set  M = d\  + d2  -I t-  dq,  then  q < M <t,  and  we 

have  q = M = t if  and  only  if  Q is  an  abelian  group. 

Let  fi  be  a bounded  domain  in  1R2  with  a piecewise  smooth  boundary  T,  invariant 
with  respect  to  Q,  i.e.,  sent  onto  itself  by  the  congruences  of  Q.  Also  the  boundary  F is 
invariant  with  respect  to  G,  i.e.,  for  any  7^  € G and  x € T,  (jex)  6 F. 
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Let  W(r)  be  the  real  vector  space  of  real  functions  defined  on  F.  We  can  associate 
to  any  element  7,;  of  Q a linear  transformation  T,  defined,  for  any  v £ W(r),  by 

(Tiv)(x)-.=  v{~f~lx)  leT,  (2.2) 

where  % is  a linear,  invertible  transformation  from  W(r)  onto  W(r),  and  Tj  is  the 
identity. 

Definition  2.1  A subset  V(r)  of  W(F)  is  said  to  be  invariant  with  respect  to  Q (or 
Q -invariant)  if  for  any  v £ V(r)  and  any  7*  £ Q,  TiV  £ V(r). 

Obviously  if  v is  a function  of  W(r),  not  identically  equal  to  zero,  the  set  of  functions 
{TiV,  i = 1, . . . , t}  is  invariant  with  respect  to  Q. 

Definition  2.2  Let  £ be  a linear  operator  in  V(r).  We  will  say  that  £ is  invariant  with 
respect  to  Q if  for  any  u £ V(r);  £T{U  = TiCu,  i = 1, . . . , t. 

Example  2.3  Let  V(r)  be  a suitable  Sobolev  space  and  ( £f)(x ) :=  /r  K.(x,  y)f(y)dTy 
an  integral  operator  defined  on  V(r),  with  kernel  IC(x,y). 

We  have:  Ti(£f)(x)  = frK.('y~1xiy)f(y)dTv;  since  7*  £ Q is  an  isometry,  the  mapping 
2/  — > 7*2/  preserves  the  differential  element  dry.  Thus 

C(Tif){x)  = J IC( x,  y)f('r~ly)dry  = K.(x,  'yiy)f(y)dTy  . 

Then  the  integral  operator  £ is  (/-invariant  if  the  kernel  IC(x,  y)  satisfies  the  condition 
K{x,  y)  = IC^x,  liy ) for  all  x,  y £ T , i = 1, . . . , t. 

Starting  from  the  group  Q,  the  system  of  representation  (2.1)  and  the  linear  trans- 
formations Ti  defined  by  (2.2),  we  can  introduce  M linear  transformations  of  W(r), 

4 = (e=l,...,q;k=l,...,de).  (2.3) 

1 i= 1 

Owing  to  the  property  of  the  representations  (2.1),  there  holds 

Pfk  = Ptk,  PikPev=  0 if  jziZp^  = Ti-  (2-4) 

e= 1 fc=i 


The  linear  transformations  P(k,  which  will  be  called  projection  operators , determine  a 
decomposition  of  any  vector  space  V(r)  c W(r)  invariant  with  respect  to  G,  into  a 
direct  sum  of  M subspaces  Wfc(r);  V«  (F)  is  the  co-domain  of  Pek,  viewed  as  a linear 
transformation  from  V(F)  onto  itself. 

If  Q is  a non-abelian  group,  it  is  useful  to  consider  in  the  space  W(r)  further  linear 
transformations  linked  to  the  system  (2.1).  Let  {^^(7;)}  be  a representation  of  Q of 
order  dt  > 2.  Let  us  consider  dj  linear  transformations,  already  introduced  in  [1],  defined 
as  follows 


4W  _ df 

kr  t 


k,  r = 1,...,  de. 


(2.5) 
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If  k = r,  then  = Ptk. 

Definition  2.4  Let  B(-,  •)  be  a bilinear  form  from  V(r)  x V(r)  on  R.  We  will  say  that 
is  Q -invariant  if  for,  any  u,  v eV(T), 

B(TiU,  Tiv)  = B(u,  v),  i = (2.6) 

Let  V(T)  be  a Hilbert  space  and  let  us  consider  the  following  problem 

find,  a £ V(r)  : B{u,v)  = T{v)  for  all  v G V(r),  (2.7) 

where  B{-,  •)  is  continuous  and  coercive,  and  P(-)  : V(r)  R a linear  continuous 
functional.  If  T and  V(r)  are  invariant  with  respect  to  Q,  and  V(r)  = ®|=1  V^fc(r) 
is  the  decomposition  of  V(r)  defined  by  the  projection  operators  (2.3)  the  following 
fundamental  result  holds. 

Theorem  2.5  //£?( •,  •)  verifies  the  condition  (2.6)  and  are  the  projection  operators 
defined  in  (2.3),  then  the  problem  (2.7)  can  be  decomposed  into  M independent  problems; 
find  uek  € Vik  (r)  such  that 

B{uek,  vtk)  = P(vek)  for  alluffc  € V^fe(r),  £ = 1, . . . , q;  k = 1, . . . , dt . (2.8) 

The  solution  of  (2.7)  can  be  recovered  as  u = ®9e=1  W&- 

The  above  result  can  be  applied,  under  the  invariance  hypothesis,  in  the  discrete  form 
to  the  symmetric  Galerkin  BEM  scheme  if  we  choose  the  finite  dimensional  subspace 
Uh,p(T)  defined  in  (1.5),  to  be  Q-invariant  too,  and  therefore  decomposable  as  Uh,p{T)  = 
®$=1  Uhfpi r)-  Then  the  symmetric  Galerkin  boundary  element  problem  can  be 

decomposed  into  M independent  problems  which  have  reduced  dimension  with  respect 
to  the  original  one  and  which  can  be  solved  on  parallel  processors.  Now  one  has  to 
construct  boundary  element  basis  functions  for  each  subspace  U^p(Tj.  With  some  simple 
geometries  (and  groups  of  congruences)  this  can  be  done  directly,  but  in  many  cases  this 
is  a difficult  task.  We  solve  it  here  by  applying  restriction  matrices , which  we  introduce 
in  the  next  sections,  to  the  basis  of  UhlP(T),  ordered  as  a vector.  Since  there  is  a ono- 
to-one  correspondence  between  the  standard  boundary  element  shape  functions  and  the 
nodes  of  the  mesh  fixed  on  F,  in  the  following  we  will  work  directly  on  the  nodes  of  the 
boundary. 

3 Elementary  restriction  matrices 

In  this  section  we  introduce  suitable  matrices  depending  only  on  the  group  Q and  on  the 
system  of  representations  (2.1),  which*will  be  called  elementary  restriction  matrices.  In 
the  following  sections  we  will  see  how,  starting  from  these,  we  can  construct  restriction 
matrices  relative  to  a mesh  defined  on  T.  We  fix  a finite  group  Q = {7!,...,  74}  of 
congruences  of  Rm  and  a system  (2.1)  of  orthogonal  irreducible,  pairwise  inequivalent 
representations  of  Q.  Q always  admits  the  representation  (1,  1, . . . , 1}  which  we  indicate 
by  (ti/1)  (7* ) } ; let  us  order  the  remaining  representations  (2.1)  with  increasing  order  df, 
let  {wl1l(7j)}, . . . , {u>ls)(7j)}  be  the  representations  of  order  1.  If  Q is  an  abelian  group 
one  has  s = q = t and  d\  = • = dt  = 1.  If  Q is  a nonabelian  group,  it  holds 

s < q < t and  therefore  d\  = d2  = ■ • ■ = ds  = 1 , 2 < ds+\  < ■■■  <dq. 
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Let  Q be  an  abelian  group.  We  will  call  elementary  restriction  matrices  the  following 
t matrices,  with  1 row  and  t columns 

Rn  = -j=(uW('h)  •••  w(f)(7t))  , £ = l,...,t.  (3.1) 


Since  representations  {lj^  (7*)}  are  real,  it  follows  that  0/^(7,;)  = ±1,  for  £,  i = 1, . . . , t. 

Let  Q be  a nonabelian  group.  Correspondingly  to  the  representations  {0/^(7*)}  of 
order  1 of  the  system  (2.1),  we  introduce  matrices  R{\  with  1 row  and  t columns 

Rn  = ~(ww(7i)  ...  ww(74))  , l = \ (3.2) 

We  obtain,  in  this  case,  s matrices.  Let  now  {^^(7^)}  be  a representation  of  the  system 
(2.1)  of  order  de,  with  de  > 2.  With  k = 1 de  fixed,  let  us  consider  the  following 
matrix,  with  d(  rows  and  t columns 


Reu  = 


l 4k  (7 1) 
i ^(7 1) 


V 4?k(7  1) 


4fc(72)  • 

4k(7t)  ^ 

4k  (72)  •• 

4k  (7t) 

4i( 72)  ■■ 

4?fc(7<)  ) 

(3.3) 


Due  to  the  orthogonality  properties  of  the  representation {u/b1  (7^ ) } j matrix  Ra  has  pair- 
wise orthonormal  rows.  Therefore  the  rank  of  matrix  Rek  is  d(.  For  any  representation 
{0/^(7*)}  we  obtain  de  matrices  R(k  ( k = 1, . . . , de).  Matrices  Reu  (£  — 1, q;  k = 
1,.. de)  defined  in  (3.2)  and  (3.3)  will  be  called  elementary  restriction  matrices.  The 

total  number  of  these  matrices  is  M,  with  M = d\  +d2-\ \-dq.  The  matrices  defined  in 

(3.1)  or  (3.2)-(3.3)  satisfy  some  properties,  easily  deducible  from  orthogonality  relations 
(2.4)  and  which  we  summarise  in  the  following. 


Theorem  3.1  ([1])  The  M elementary  restriction  matrices  defined  by  (3.1)  or  (3.2)- 
(3.3)  verify  the  relations 

<1  dc 

RtkRlk  = Ide , RekR},k,  = 0 if  (£,  k)  £ (£',  k') , EE  Rh-Rek  = I (3-4) 

e=i  fc= 1 

where  Idl , 1 are  identity  matrices  of  order  de  and  t respectively. 


4 'H(Sa)  spaces  and  elementary  restriction  matrices 


Let  T be  the  piecewise  smooth  boundary  of  ft,  invariant  with  respect  to  Q.  and  a £ F. 
Consider  the  ordered  set 

Da  = {0,72  1a,...,7fi"1o},  (4.1) 

and  the  space  7f(Ea)  of  real  functions  defined  in  En.  A natural  basis  B in  W(Ea)  is 
formed  by  functions  having  value  1 in  a point  of  Ea  and  0 in  the  remaining  points. 
Having  indicated  with  x the  function  of  B with  value  1 in  the  point  a,  we  obtain  the 
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ordered  basis  B = {x(x)-,  x{l2%),  ■ • • > xiltx)},  such  that,  of  course, 

U(T,a)  = span{x(a;),x(72x), . . • ,x(7 tx)}.  (4.2) 

H(Ea)  is  a vector  space  with  finite  dimension  n < t,  invariant  with  respect  to  Q (since 
£0  is  invariant  with  respect  to  Q)  and  therefore  decomposable  into  direct  sum  of  M 
subspaces  W«c(£a).  Having  set  n g = dim  Utkina),  we  have  n = dgng. 

Definition  4.1  We  say  that  a is  a generic  point  ofT  ( with  respect  to  the  group  Q)  if 
dimTf  (£a)  = t or,  equivalently,  if  all  the  elements  of  £a  are  distinct. 

The  following  results  hold. 

Theorem  4.2  ([1])  Having  fixed  any  point  a G T,  if  {wW(7j)}  is  a representation  of 
order  1,  then  Hgi(iua)  = span{Pax}  and  ng  < 1.  // {0/^(7;)}  is  a representation  of 
order  dg  > 2,  one  has 

Heki^a)  = span{^x, • • • , k = l,...,de,  (4.3) 

and  therefore  ng  < dg.  If  a is  a generic  point,  then  ng  — dg  for  any  i. 

Let  now  V 4 be  the  column  vector  (x(x)"X{X‘2x),  ■ ■ ■ , x(ltx))* , whose  order  is  related  to 
that  one  fixed  for  the  elements  of  Q.  Corresponding  to  the  representations  of  order  1 of  Q, 
for  the  elementary  restriction  matrices  defined  in  (3.1),  (3.2)  we  have  RgiV1  = \ftPg\X- 
From  Theorem  4.2,  it  follows  that 

Hn(Xa)=spm{Rg1Vt}.  (4.4) 

Corresponding  to  the  representations  of  order  dg  > 2,  for  the  elementary  restriction 
matrices  defined  in  (3.3)  we  have  RgkV 1 — y/ 1 / dg^A^lx,  A ■ From 

(4.3),  it  follows  that 

Hgk(£a)  = span {Rg^}  . (4.5) 

In  both  cases,  if  a is  a generic  point,  the  components  of  the  vector  RgkV1  constitute  a 
basis  in  HgkifBa)-  Therefore,  for  any  generic  point  a,  the  elementary  restriction  matrix 
Rgk  represents  the  projection  operator  Pgk  from  7Y(£0)  onto  Hgk^a),  if  we  choose  V1 
as  a basis  in  7f(£0)- 

Now,  we  want  to  construct  elementary  restriction  matrices  Rgk  which  represent  the 
projection  operators  Pgk  from  'HifRa)  onto  Hgki'Ba)  for  nongeneric  points.  Therefore  let 
us  suppose  a to  be  a nongeneric  point,  i.e.,  such  that  the  functions 

X(z),x(72z),--.,x(7tz)  (4-6) 

are  linearly  dependent.  Let  n be  the  maximum  number  of  linearly  independent  functions 
among  (4.6)  and  let  the  following  functions  be  linearly  independent, 

X{lhx),---,x{linx)-  (4-7) 

It  is  convenient  to  order  the  functions  (4.7)  with  increasing  index  ia;  therefore  let  us 
suppose  i\  < *2  < • • • < in-  In  this  case  elementary  restriction  matrices  Rgk  will  have  n 
columns.  The  number  ng  of  rows  (ng  < dg)  of  each  Rgk  is  not  determined  by  ii,  V,  - ■ ■ , in- 
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In  general,  we  only  can  say  that  matrices  Rp\ , . . . , Rf  d,  have  the  same  number  rip  of  rows, 
where  m = dim  Utkina)- 

Then  we  now  consider  a significant  class  of  nongeneric  points.  Having  fixed  £(£  = 
2, . . . , t),  let  Ir(r)  be  the  set  of  all  points  a £ T such  that 

a = (4.8) 

From  (4.8)  it  follows,  for  any  i : x(lix)  = This  implies  that  the  functions  (4.6) 

are  naturally  subdivided  into  subsets  and  any  subset  contains  coincident  functions.  Then 
we  can  obtain  elementary  restriction  matrices  for  the  space  Ti (£a)  with  a £ /f  (T)  starting 
from  elementary  restriction  matrices  built  in  Section  3,  with  the  following  procedure, 

• Let  us  sum  to  each  column  of  index  ia  (a  = 1, . . . , n)  all  the  columns  of  index  j, 
with  j such  that  7 ~la  = 7 “xa.  We  indicate  with  Rpk  the  obtained  matrices,  all 
with  de  rows  and  n columns,  but  not  all  full-rank  matrices;  some  of  these  may  be 
zero  matrices. 

• Let  us  extract  from  nonzero  matrices  Rpu  submatrices  Rpk  made  up  of  np  linearly 
independent  rows. 

• Finally,  let  us  construct  from  Rpp  matrices  Rpk  with  a row-orthonormalization  pro- 
cedure. 

The  (nonzero)  matrices  Rpk  verify  the  properties  expressed  by  Theorem  3.1.  Further- 
more, matrices  Rek,  applied  to  the  vector  Vn  = (x(^ixx), . . . , x(7inx))*  corresponding 
to  a point  a £ Ie(T).  give  vectors  whose  components  constitute  a basis  for  7 i«-(£a)- 
For  this  reason  they  represent  the  projection  operators  from  7-f(£0)  onto  7f^(£a),  for 
any  a £ J*(r).  Then  we  will  say  that  the  matrices  Rpk,  with  np  rows  and  n columns, 
are  elementary  restriction  matrices  for  the  space  W(£a)  relative  to  points  a £ Ie(T). 
Furthermore  n = dene. 

5 'W(S)  spaces  and  restriction  matrices 

Let  r be  the  piecewise  smooth  boundary  of  Q,  £ a set  formed  by  N points  of  T consti- 
tuting a not  necessarily  uniform  mesh  defined  on  F.  Let  us  suppose  F and  £ invariant 
with  respect  to  Q.  Let  W(£)  be  the  vector  space  of  real  functions  defined  in  £.  7f(£)  is 
a iV-dimensional  vector  space,  invariant  with  respect  to  Q\  this  is  due  to  the  fact  that  £ 
is  invariant  with  respect  to  Q.  A natural  basis  B in  7f(£),  invariant  with  respect  to  Q, 
is  formed  by  functions  having  value  1 in  a point  of  £ and  0 in  the  remaining  points.  In 
order  to  more  easily  construct  restriction  matrices  for  the  space  7T  ( £) , or  equivalently 
for  the  mesh  £,  it  is  suitable  to  introduce  in  the  set  £ the  following  equivalence  relation. 

Definition  5.1  We  say  that  a point  a1  is  equivalent  to  a"  if  there  exists  an  element 
7 i £ G such  that  a"  = 7 ~xa!  (and  therefore  a'  = 7,0"). 

The  points  of  the  set  £ are  then  subdivided  into  r equivalence  classes.  If  r = 1 one 
has  7f(£)  = 7i(£a),  with  a £ £.  Then  let  us  suppose  r > 2.  We  order  the  points  of  the 
set  £ as  follows;  having  indicated  with  a\, . . . ,ar  r pairwise  inequivalent  points  of  £,  we 
consider  the  following  ordered  points 

ai,72  lai.---i7r1«i.«2,7j102,---,7rla2,---,ar,7jVI---)7r10r-  (5-1) 


Geometrical  symmetry 


If  points  (5.1)  are  distinct,  we  have  N = rt.  If  some  points  among  (5.1)  coincide,  we  will 
erase  from  the  sequence  (5.1)  a point  if  it  is  equal  to  a previous  one.  Then  a sequence 
of  N points,  with  N < rt , will  remain,  with  points  equivalent  to  cq,  n<2>  equivalent 
to  a2,. . . , n equivalent  to  ar.  In  both  cases  7i(E)  = W(SQl)  ©W(Sa2)  ® • ■ • ® ft(Ear) , 
with  dim7f(£aj.)  = < t,  j = 1, . . . ,r  and  N = n ^ + n ® + • • • + n^r\  We  indicate 

by  C$  the  elementary  restriction  matrices  relative  to  the  space  constructed  as 

indicated  in  Section  4.  Let  np  be  the  number  of  rows  of  the  matrix  C'P ; having  fixed  j, 
the  number  of  columns  of  matrices  Cj%,  for  any  £ and  k,  is  n^\  We  consider  therefore 
the  following  M block  matrices 


fdS 


Rtk  = 
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O \ 
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(5.2) 


\ 0 o o ...  cil>  J 


with  N£  = nP  +nP  H h rSp  rows  and  N columns,  from  which  we  have  to  eliminate  the  ■ 

possible  zero  rows.  Matrices  Rik  determined  by  this  procedure,  which  we  call  restriction 
matrices  for  the  space  7f(S)  of  dimension  N,  have  rank  equal  to  the  number  Ne  of  the 
remained  rows  and  for  these  matrices  properties  expressed  in  Theorem  3.1  still  hold.  In 
both  cases,  we  have  the  following  theorem. 


Theorem  5.2  Considering  the  basis  B in  W(E)  as  a column  vector  VN  with  the  or- 
der deduced  from  (5.1),  the  components  of  the  vector  RekVN  form  a basis  in 
Therefore  the  M matrices  R ek,  having  fixed  in  W(£)  the  ordered  basis  VN , determine  a 
decomposition  of  H (E)  in  M subspaces,  which  coincides  with  the  one  obtained  with  the 
projection  operators  P ek 


Preliminary  numerical  results  appear  promising;  algorithms  for  potential  and  linear 
elasticity  problems  are  being  implemented  on  parallel  processors  to  analyse  the  efficiency 
of  the  proposed  approach. 
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